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^ ; ABSTRACT 

H . We report results from a 1 week multi-wavelength campaign to monitor the BL Lac object 

S5 0716+714 (on December 9-16, 2009). Nine ground-based telescopes at widely separated 
longitudes and one space-based telescope aboard the Swift satellite collected optical data. 
Radio data were obtained from the Effelsberg and Urumqi observatories, X-ray data from 
Swift. In the radio bands the source shows rapid (0.5 — 1.5) day) intra-day variability 
with peak amplitudes of up to ^ 10%. The variability at 2.8 cm leads by about 1 day the 
variability at 6 cm and 1 1 cm. This time lag and more rapid variations suggests an intrinsic 
contribution to the source's intraday variability at 2.8 cm, while at 6 cm and 1 1 cm interstellar 
scintillation (ISS) seems to predominate. Large and quasi-sinusoidal variations of ^ 0.8 mag 
were detected in the V, R and I-bands. The X-ray data (0.2-10 keV) do not reveal significant 
variability on a 4 day time scale, favoring reprocessed inverse-Compton over synchrotron 
radiation in this band. The characteristic variability time scales in radio and optical bands 
are similar. A quasi-periodic variation (QPO) of 0.9 — 1.1 days in the optical data may be 
present, but if so it is marginal and limited to 2.2 cycles. Cross-correlations between radio and 
optical are discussed. The lack of a strong radio-optical correlation indicates different physical 
causes of variability (ISS at long radio wavelengths, source intrinsic origin in the optical), 
and is consistent with a high jet opacity and a compact synchrotron component peaking at 
~ 100 GHz in an ongoing very prominent flux density outburst. For the campaign period, 
we construct a quasi-simultaneous spectral energy distribution (SED), including 7— ray data 
from the FERMI satellite. We obtain lower limits for the relativistic Doppler-boosting of 5 ^ 
12 — 26, which for a BL Lac type object, is remarkably high. 
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1 INTRODUCTION 

Blazars are the extreme subset of active galactic nuclei that are 
usually taken to include both BL Lacertae objects (BL Lacs) and 
flat spectrum radio quasars (FSRQs) because both are character- 
ized by rapid variability across the electromagnetic spectrum and 
significant radio to optical polarization. Their spectra are gener- 
ally well modelled by synchrotron and inverse Compton emission 
from relativistic jets aligned nearly (< 10°) with the line of sight 
(e.g., Urry & Padovani 1995). Blazar flux variations are seen on 
timescales ranging from a few minutes through days and months 
to decades. Blazar variability timescales can be divided into three 
classes: changes from minutes to less than a day are variously 
called microvariability, intra-night variability, or intra-day variabil- 
ity (IDV), which is the term we shall use; those from a few days to a 
few weeks are usually known as short timescale variability (STV); 
flux changes from months to many years are called long term vari- 
abiUty (LTV; e.g., Gupta et al. 2004). 

The vast majority of optical observations of blazars take place 
over single nights or several nights in succession at a given tele- 
scope, but these cannot provide the continuous coverage one would 
like to have in order to properly characterize and understand IDV 
and STV. The possibility that some coherent, even if temporary, 
fluctuations are present in blazar light curves is very important 
to investigate. Certainly one improves the chance of seeing such 
variations by having a long, densely sampled, light curve. These 
requirements have led to several intensive campaigns in the past 
that have looked at individual blazars with multiple telescopes 
at widely separated longitudes over a few days through several 
months where simultaneous observations are obtained at different 
bands of the electromagnetic spectrum (e.g., Villata et al. 2000, 
2006, 2008, 2009a, 2009b; Raiteri et al. 2003, 2005, 2007, 2008a, 
2008b, 2008c; Bottcher et al. 2005, 2009; Ostorero et al. 2006; 
Agudo et al. 2006; Fuhrmann et al. 2008; Larionov et al. 2008, 
and references therein). 

We report results of a 7 day long multi-telescope campaign on 
S5 0716+714 with a focus on IDV in the optical and at some radio 
bands. The campaign was conducted on December 9 - 16, 2009 
with the most dense sampling and best frequency coverage during 
December 11-15. There have been many studies of this blazar, 
mainly in optical bands, since it is bright, at a high declination and 
apparently always quite variable in the visible (e.g., Quirrenbach 
et al. 1991; Wagner et al. 1996; Sagar et al. 1999; Wu et al. 2005, 
2007; Montagni et al. 2006; Stalin et al. 2006; Pollock et al. 2007; 
Gupta et al. 2008; Poon et al. 2009; Rani et al. 2010a, 201 1 ; Chan- 
dra et al. 2011, and references therein). Our observing campaign 
also inckidcs flux density measurements made in the centimeter 
and millimeter radio, optical and X-ray bands. 

Theoretical models that seek to explain optical intra-day vari- 
ability in AGN invoke several different mechanisms involving the 
accretion disk, including pulsation of the gravitational modes of 
the gaseous disk (e.g., Kato & Fukue 1980; Nowak & Wagoner 
1992) or orbital signatures from "hot-spots" in the gas surround- 
ing the black hole, either from the disk itself or the corona above 
it (e.g., Zhang & Bao 1991; Mangalam & Wiita 1993). However, 
for blazars, particularly in high states, the variability almost cer- 
tainly arises within the Doppler boosted relativistic jets and may 
well result from relativistic shocks in the jet (e.g. Marscher & Gear 
1985) or turbulence behind such shocks (e.g. Marscher, Gear & 
Travis 1992; Marscher et al. 2008), from helical motion (e.g. Qian 
et al. 1991, Camenzind & Krockenberger 1992), instabiUties (e.g. 



Hardee et al. 2005), or slight changes in viewing angles (e.g. Gopal- 
Krishna & Wiita 1992). 

A key motivation of this campaign was to look for possible 
correlations between radio and optical variability. A positive cor- 
relation could substantially constrain the origin of the variability 
in the radio bands and would strongly favor an intrinsic origin 
over one due to interstellar scintillation (e.g. Simonetti, Cordes & 
Heeschen 1985, Rickett 1990). There have been earlier papers in 
which the variability properties of S5 0716-1-714 were discussed 
in both of these bands (e.g., Quirrenbach et al. 1991; Wagner & 
Witzel 1995; Wagner et al. 1996; Ostorero et al. 2006; Fuhrmann 
et al. 2008). Quirrenbach et al. (1991) reported correlated opti- 
cal and radio variability in the source which, however, was not 
again seen in later observational campaigns. Optical variability was 
claimed to be associated with changes in the radio spectral index 
(Qian et al. 1995 & 1996). This source is known to vary on dif- 
ferent timescales (IDV to STV to LTV) and has been claimed to 
exhibit quasi-periodic oscillations (QPOs) on all these timescales. 
On IDV timescales 0716+714 has shown QPOs on various occa- 
sions with timescales ranging from ~ 15 min to 73 min (Gupta et 
al. 2009; Rani et al. 2010b). On STV timescales there is weak evi- 
dence for quasi-periodicity appearing on timescales of ~ 1 day and 
~ 7 days (Quirrenbach et al. 1991; Wagner 1992; Heidt & Wag- 
ner 1996). On LTV timescales claims have been made for possible 
near-periodicities of ^ 3.0+0.3 years (Gupta et al. 2008). We adopt 
a redshift of ^; = 0.31 for this source (Nilsson et al. 2008) but the 
lack of spectral lines means that it is still quite uncertain (±0.08). 

In §2 we describe the observations and data reductions. We 
present our results in §3 and §4 includes a discussion and our con- 
clusions. 



2 OBSERVATIONS AND DATA REDUCTIONS 
2.1 Radio data 

The source S5 0716-1-714 was observed on December 11-15, 2009 
with the MPIfR 100m telescope at Effelsberg, Germany and with 
the 25m Nanshan radio telescope of the Urumqi Observatory, 
China. The observations and source selection at both observato- 
ries were coordinated and done in a similar manner. At Effelsberg 
the observations were performed at 3 frequencies using the sec- 
ondary focus heterodjTie receivers operating at 2.7, 4.85 and 10.5 
GHz. At Urumqi the observations were done using a Cassegrain 
focus receiver, which operates at a frequency of 4.8 GHz. At both 
telescopes the flux density measurements were made in a similar 
way, with repeated cross-scans in azimuth and elevation. Frequent 
switching between target and calibrator sources on time scales of a 
few minutes allowed us to monitor the gain variations introduced 
by the receiving system and the atmosphere. The gain variations 
seen in the calibrator sources were used to improve the flux den- 
sity calibration and correct for time and elevation dependencies. 
This observing and calibration technique is well established and 
has been applied before for various intra-day variability observa- 
tions. The details of the observing strategy, the receiver parameters 
and the calibration methods have been described before (Gabanyi 
et al. 2007; Fuhrmann et al. 2008; Marchili et al. 2010). As sec- 
ondary calibrators the two non- variable radio sources 0951-1-699 
and 0836-1-710 were used. The absolute flux density scale was set 
using 3C 286 and NGC 7027, adopting the scale of Zijlstra et al. 
(2008). 
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2.2 Optical data 

Our official campaign on S5 0716+714 began on 11 December 
2009 and ran through 15 December 2009. Additional data was ob- 
tained on 9, 10, 16 and 17 December, though the coverage was not 
as dense. We briefly describe the telescopes and cameras that were 
involved in these observations. We note that as this source is point- 
like the different apertures used at different telescopes have negli- 
gible effects on the measured fluxes and agreement was excellent 
whenever different telescopes provided overlapping data. 

The source was observed on 12 December 2009 with the 23.5 
cm f/10 Schmidt-Cassegrain telescope located on the roof of the 
Department of Physics of La Sapienza University (Roma, Italy), 
with a Bessel R filter and a CCD camera with a Peltier cooled Ko- 
dak KAF 1603ME CCD chip. Due to the lack of a guiding sys- 
tem, images were taken in sequences of 9-10 frames with 60 s 
exposures and then stacked and summed together to increase the 
S/N ratio. Sky flats were used for flat fielding. Data reduction and 
analysis were performed with standard Image Reduction and Anal- 
ysis Facility ( I RAF )[3 routines. Aperture photometry was made 
with a 3.3 arcsec (2 pixels) radius. Simultaneous observations were 
made with the 31 cm f/4.5 Newtonian telescope located at Greve in 
Chianti (near Florence, Italy) equipped with a CCD camera with a 
back-illuminated Peltier cooled Site SIA502A chip and B, V, R, I 
Bessel filters. A four color sequence was made at the beginning of 
observations and then R filter monitoring began with 180 s expo- 
sure times. I RAF standard tasks were used for data reduction and 
analysis; aperture photometry was made with a 7 arcsec (2 pixels) 
radius. Photometric errors were estimated from the rms deviations 
of the reference stars, and found to be generally less than 0.01 mag. 
The light curves of the two instruments were in very good agree- 
ment. Reference standard stars were taken from Villata et al. (1998) 
in the same blazar field. 

At the MDM Observatory on the south-west ridge of Kitt 
Peak, Arizona, USA, data were taken for limited periods during 
the nights of 9, 10, 11 and 12 December with the 1.3 m McGraw- 
Hill Telescope, using the Templeton CCD with B, V, R, and I filters. 
The standard data reduction, using I RAF, included bias subtraction 
and flat-field division. Instrumental magnitudes of S5 0716+714 
plus four comparison stars in the field (Villata et al. 1998) were 
extracted using the IRAF package DAOPHOT0 (Stetson 1987) with 
an aperture radius of 6 arcsec and a sky annulus between 7.5 and 
10 arcsec. 

The observations at the Abastumani Observatory were con- 
ducted on 9, II, 12, 14, 15 and 16 December 2009 at the 70-cm 
meniscus telescope (f/3). These measurements were made with an 
Apogee CCD camera Ap6E (IK x IK, 24 micron square pixels) 
through a Cousins R filter with exposures of 60-120 sec. Reduction 
of the image frames were done using DAOPHOT II. An aperture 
radius of 5 arcsec was used for data analysis. 

Observations of S5 0716+714 were carried out on 13 Decem- 
ber 2009 using the 50/70-cm Schmidt telescope at Rozhen National 
Astronomical Observatory, Bulgaria. The telescope is equipped 
with a FLI (Finger Lakes Instrumentation) ProLine PL 1 6803 CCD 
(4096 pixels x 4096 pixels), BV Johnson and RI Cousins filters. 
Standard data reduction used MIDAS and included bias subtraction 



IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 

^ Dominion Astrophysical Observatory Photometry software 



and flat-field division. Instrumental magnitudes of S5 0716-1-714 
and comparison stars in the field (Villata et al. 1998) were extracted 
using the MIDAS package DAOPHOT with an aperture radius of 5.6 
arcsec (2xFWHM). 

Observations of the source on 11, 13 and 15 December 2009 
were carried out using the 1.04 m Sampumanand telescope (ST) 
located at Nainital, India. It has Ritchey-Chretien optics with a 
f/13 beam equipped with a CCD detector which is a cryogenically 
cooled 2048 pixel x 2048 pixel chip mounted at the Cassegrain 
focus and Johnson UBV and Cousins RI filters. Each pixel of 
the CCD chip has a dimension of 24 firn^, corresponding to 
0.37arcsec^ on the sky, thereby covering a total field of ~ 13' x 
13'. Data processing (bias correction, flat-fielding and cosinic ray 
removal) was done by the standard routines in IRAF and pho- 
tometry of the blazar and standard stars in its field employed a 
stand-alone version of DAOPHOT II. Aperture photometry was 
carried out with four concentric aperture radii, i.e., ~ IxFWHM, 
2xFWHM, 3xFWHM and 4xFWHM. We found that aperture 
radii of 2xFWHM almost always provided the best S/N, so we 
adopted that aperture for our final results. 

Observations also were carried out on 13 December 2009 us- 
ing the 61-cm Boiler and Chivens reflector at Sobaeksan Opti- 
cal Astronomy Observatory in Korea. The FLI CCD camera with 
thermo-electric cooled Kodak KAF 4301 2K CCD chip set and 
standard R Cousins filter were used for the observations. All im- 
ages were reduced using standard IRAF tools. Aperture photome- 
try parameters that maximized the S/N were an aperture radius of 
1.5 X FWHM, an inner radius of the sky annulus of 5 x FWHM, 
and a sky annulus width of 10 pixels. 

On 12, 13 and 14 December 2009 the source was observed at 
the 1.2 m Cassegrain telescope at the Michael Adrian Observato- 
rium of Astronomic Stiftung Trebur, Germany. A Roper Scientific 
EEV 1340 EB CCD camera and a Cousins R filter were employed. 
The data were reduced using MIRA Pro 7 software. An aperture 
radius of 2.9 arcsec, and inner and outer sky annulus radii of 10.4 
and 13.8 arcsec, respectively, were used. 

Observations of S5 0716+714 were carried out on 12 Decem- 
ber 2009 at the 0.9 m optical SARA (Southeastern Association for 
Research in Astronomy) telescope at Kitt Peak National Obser- 
vatory, USA. The camera was an Apogee U42 with Johnson and 
Cousins UBVRI filter set. MIRA software was used for image pro- 
cessing and data analysis. The data were analyzed with aperture 
size of 4 arcsec. The blazar light curve was calibrated using the 
local standard stars present in the blazar field (Villata et al. 1998). 

The final set of ground based optical observations of 
S5 0716+714 were carried out on 10 and 16 December 2009 at 
the 2.56 m Nordic Optical Telescope (NOT), Canary Islands, Spain 
using ALOFSC (10 Dec) and MOSCA (16 Dec) using UBVRI fil- 
ters. The data were reduced using standard IRAF procedures, in- 
cluding those for debiasing and twilight flat-correction. The pho- 
tometry was done using IRAF/APPHOT with the aperture radius 
chosen to be close to the FWHM of the image. 

S5 0716+714 also was observed by the Swift satellite's 
Ultraviolet-Optical Telescope (UVOT, Roming et al. 2005) in the 
V-band during 54 pointings conducted between December 1 1 and 
15, 2009. The pointings were typically separated by the 96 minutes 
period of Swift's orbital revolution. Observations from the space- 
based platform have the obvious advantage of not being interrupted 
by the day/night cycle or unfavorable weather, and therefore can 
probe variability on timescales of hours to days with dense and uni- 
form sampling. The UVOT is a D = 300 mm, f = 38 10 mm modified 
Ritchey-Chretien telescope equipped with a micro-channel plate in- 



4 Alok C. Gupta et al. 



tensified CCD detector operated in photon counting mode. Detec- 
tors of tliis type may provide information about the time of arrival 
of individual photons and they have no specific saturation limit, 
but these advantages come at the cost of a nonlinear response to 
the number of incoming photons because of the dead time after 
each registered photon event known as the coincidence loss (or 
pile-up) effect (Poole et al. 2008, Breeveld et al. 2010). The VaST 
software (Sokolovsky and Lebedev 2005) based on the SExtrac- 
tor code (Berlin and Arnouts 1996) has been applied to conduct 
aperture photometry of the UVOT images. The VaST software is 
designed to deal with imaging data obtained with non-linear detec- 
tors (e.g. Kolesnikova et al. 2008, 2010) and has been successfully 
applied for the UVOT data reduction before (Sokolovsky 2009). 
The magnitude scale was set using comparison stars 4-8 from Vil- 
lata et al. (1998). The data reduction technique we employed al- 
lows us to avoid uncertainties in the coincidence loss correction for 
bright sources and enables direct comparison of the UVOT results 
with those obtained with ground-based telescopes after a R-V color 
correction of 0.43 ± 0.04 mag is applied (Rani et al. 2010a). 

2.3 X-ray data 

Swift's X-ray telescope (XRT, Burrows et al. 2005) observed the 
source simultaneously with the UVOT in the 0.2-10 keV energy 
range. The XRT was operating in the photon counting (pc) mode. 
The total exposure of 15.5 ksec was collected during 54 satellite 
pointings performed between 2009 December 11-15; each point- 
ing was about 300 sec long. Observations were processed using the 
xrtpipeline tool from the HEASOFT v6 . 8 package applying 
the standard filtering criteria. The mean count rate during the ob- 
servations was 0.53 ± 0.01 cts/sec which allows us to neglect the 
pile-up effect. 

For the spectral analysis, all sub-exposures were combined 
into a single event file using XSELECT, the corresponding expo- 
sure maps were generated with xrtexpomap and combined us- 
ing XI MAGE. The source spectrum was extracted from circular 
region of 20 pixel radius. The background counts were extracted 
from a region away from the source. The auxiliary response file 
(ARF) was generated with xrtmkarf, and the response matrix 
swxpcOtol2s 6_2 007 90 IvOl 1 . rmf was used. The spec- 
trum was re-binned with the tool grppha to contain at least 25 
counts per energy bin to enable use of the statistic. It was mod- 
eled with an absorbed power law, taking the neutral hydrogen col- 
umn density fixed to the Galactic value of A'hi = 2.0 x 10^" cm~^ 
from Kalberla et al. (2005) while the photon index, F, and the nor- 
malization were left as free parameters. The spectral modeling was 
conducted in XS PEC vl2.5.1n. 

For the light curve analysis, the same procedure was applied 
to each individual sub-exposure. The value of F was fixed to that 
obtained from the combined spectrum to decrease the uncertainty 
of the flux measurements. Unabsorbed fluxes were computed by 
integrating the power law model in the 0.2-10 keV energy range. 



3 ANALYSIS AND RESULTS 

3.1 Variability curves: general trends and behavior 

3.1.1 Radio 

We characterize the variability parameters in the same way as in 
Fuhrmann et al. (2008), by the variability index, m, the noise-bias 
corrected variability amplitude, Y, and a reduced chi-square value. 



Xr, for a fit to a constant flu}|f[ In Table 1 we summarize the re- 
sults. In column I we give the observing frequency, in column 2 the 
mean flux density and its error. In column 3 the variability index of 
0716+714 is given, while the next column contains the variabil- 
ity index of the secondary calibrators (0951-1-699 and 0836-1-710), 
which is a measure of the residual calibration errors. Column 5 has 
the variability amplitude Y, followed by the reduced Xr in column 
6, the number of measurements in column 7, the reduced Xr value 
corresponding to a significance level of 99.9 per cent in column 8 
and the observing telescope in the last column. 

As it is seen from Fig.[T] the overall variability is mild, with a 
decreasing trend of 5S/ < S > of ~ 8 per cent at 2.7 and 4.8 GHz, 
and a ~ 7 per cent increase in SS/ < S > toward the end of the ex- 
periment at 10.5 GHz. At all frequencies the source varies on time 
scales of 0.5-1 days with a variability index m of about 1-2 per 
cent. Since the influence of the atmosphere and residual calibration 
errors increase with frequency, we also see a slight increase of the 
variability index of the secondary calibrators (mo) with increas- 
ing frequency, from 0.25-0.33 per cent at 2.7 GHz and 4.8 GHz to 
l.I per cent at 10.5 GHz. The calibration-bias corrected variabil- 
ity index Y takes into account the calibration uncertainties in mo 
and their frequency dependence. We therefore see a decrease of the 
variability amplitude Y for 0716-1-714 from 7 per cent at 2.7 GHz 
to 4.9 per cent at 10.5 GHz. Thus the rms-amplitude of the intra- 
day variability of 0716+714 decreases with increasing frequency, 
as it could be expected for (weak) refractive interstellar scintillation 
(e.g. Walker 1998, Beckert et al. 2002). We however note that our 
present finding is opposite to an observed increase of the variabil- 
ity amplitude with frequency for inter-day variability in this source 
when it was observed in November 2003 (Fuhrmann et al. 2008). 
Temporal changes of the frequency dependence of the RADIO vari- 
ability index in 0716-)-714 are well known and can be interpreted as 
a result of variable source size and the related change of the relative 
dominance of interstellar scintillation over source intrinsic variabil- 
ity (cf. Krichbaum et al. 2002). 

We qualitatively compare the observed modulation indices 
with those obtained from an analytical solution for the strength 
of interstellar scintillation (ISS) in the weak, quenched scattering 
regime (using equations of Beckert et al. 2002; see also Fuhrmann 
et al. 2008). For a plasma screen of 1 pc thickness and a turbulence 
strength similar to the value obtained from pulsar scintillation in 
the Local Bubble, the observed variability indices shown in Table 
1 can be reproduced by assuming a scintillating component of size 
6* ~ 0.1 - 0.2 mas (VLBI core size at 5 GHz, Bach et al. 2006) and 
an adopted screen distances of ~ 50 — 200 pc for the Local Bubble 
(Bhat et al. 1998). We further assume that the scintillating source 
component contains about 70 per cent of the source's total flux den- 
sity as indicated from the VLBI compactness and re-scale the ob- 
served modulation indices accordingly. Figure |2] shows the scaled 
and noise bias corrected variability index (~ \/ (m/O.?)^ — mj^) 
plotted versus frequency together with models for screens at 100 pc 
and 200 pc distances. The data are in good agreement with this ISS 
slab model, which also reproduces the observed variability time 
scale of the order of ~ 1 day (see following sections). At 10.5 GHz, 
however, the observed variability index is larger than expected from 
the model. At 4.8 GHz an annual modulation of the intra-day vari- 
ability time scale has been found which is caused by ISS (Liu et 
al. 2012). Underlying source intrinsic variability may explain this 



The .source is considered to be variable if the x^-test gives a probability 
of 0.001 for the assumption of constant flux density 



Multiwavelength Intraday Variability of the BL Lac S5 0716+714 5 



Table 1. The variability analysis parameters at radio wavelengths 



Freq. 


Mean Flux 


m 


mo 


Y 


X?. 


N 


99.9% 


Telescope 


GHz (cm) 


(Jy) 


(%) 


(%) 


(%) 










2.7 (11.0) 


1.03±0.024 


2.35 


0.30 


7.0 


63.595 


127 


1.44 


Elfelsberg 


4.85 (6.0) 


1.46±0.024 


1.67 


0.25 


4.94 


24.360 


132 


1.44 


Elfelsberg 


4.85 (6.0) 


1.46±0.029 


2.01 


0.33 


5.94 


16.256 


117 


1.45 


Urumqi 


10.5 (2.8) 


2.45±0.048 


1.96 


1.10 


4.87 


4.865 


101 


1.49 


Effelsberg 



m=variability index = as/<S>,crs standard deviation 
mo = variability index of the secondary calibrators 



Y= m^ — rriQ = bias corrected variability amplitude (see Fuhrmann et al. 2008) 

reduced Chi-square 
N= number of data points 
99 9% ~ reduced Chi-square con'esponding to a significance level of 99.9% 



excess. If we use only the first 2 days of the 2.8 cm data and by this 
omit the flux density rise seen at the end of the observations (see 
Fig. [T]l, the resulting lower variability index would agree with the 
expected frequency dependence (open and filled symbols in Fig.|2ll- 

3.1.2 Optical 

The variability is even more pronounced at optical frequencies. The 
optical light curves of the source in the B, V, R and I passbands 
shown in Fig. [TJb) display simultaneous flaring trends. Instead of 
a continuous rise or decay, multiple peaks and troughs of the flux 
were observed that give the appearance that the source might be 
showing nearly periodic variations. Although the B and I band light 
curves are less densely sampled compared to V and R, this behav- 
ior is consistent through all the optical passbands. The variability 
timescales are also comparable at all four optical passbands. De- 
tailed discussions of the variability timescales and a search for pos- 
sible quasi-periodic oscillations (QPOs) in the R-band light curve 
are given in Sections 3.2 and 3.6, respectively. 

3.1.3 X-ray 

While the source was highly active at optical frequencies no sig- 
nificant variation of X-ray flux was detected during the campaign 
(see Fig. 1(a)). The Xr test returns a 0.75 probability that the en- 
tire observed light curve is the result of random noise. From the 
visual inspection of the light curve it seems that in the beginning of 
the campaign (before JD2455 178.6) the source exhibits some vari- 
ation but it remains quiet during the second half of the campaign. 
The Kolmogorov-Smirnov (K-S) test gives only a 0.016 probabil- 
ity that flux measurements obtained before JD=2455 178.6 and af- 
ter this date are drawn from the same parent distribution. The Xr 
test returns a 0.58 probability that values before JD2455 178.6 are 
just random noise while this rises to a 0.85 probability for the val- 
ues after JD=2455 178.6. It is possible that the object exhibits small 
intraday variability in X-rays, but we were not able to detect it sig- 
nificantly with these observations; the error bars are just too large. 

3.2 Variability time scales 

To estimate any time scales present in the variability at optical and 
radio wavelengths we first employed the structure function (SF) 
analysis method. We followed Rani et al. (2009) to calculate the 
structure functions of the observed radio and optical light curves. 
Figure [3] shows the SF for the three radio light curves at 11, 6 



and 2.8 cm wavelengths. At 11 cm and 6 cm the amplitude and 
time scale of the variability are very similar, showing a break in 
the structure function (indicating a characteristic variability time 
scale) at tmin ~ 1-5 ± 0.1 day at 11 cm and tmin = 1-1 ± 0.1 
day at 6 cm. The latter time scale is visible in both, the 6 cm data 
sets from Effelsberg and Urumqi. At 1 1 cm and 6 cm another break 
near tmax = 2.9±0.1 days is only marginally significant, since it is 
longer than half of the duration of the observations. For this slower 
variability (t > 2.5 — 3 days) the data indicate higher variability 
amplitudes at 1 1 cm than at 6 cm. The SF curves at 1 1 cm and 6 
cm wavelengths are characterized by two different slopes. The SF 
curve at 11 cm follows a power-law slope, /3 — 0.90 ± 0.04 until 
tmin (with SF oc t^) and then /3 = 2.02±0.11 fromtmin totmax, 
while at 6 cm /3 = 0.98 ± 0.05 until tm^n and /3 = 1.67 ± 0.06 

from train tO tniax- 

The SF of the data at 2.8 cm, however, shows a more pro- 
nounced and faster variability than for the two longer wavelengths. 
This is also seen in the auto-correlation functions (see Section 3.3 
and Fig.|6]l. The characteristic time scales of variability at 2.8 cm 
wavelength are tmin = 0.50 ± 0.05 days and tmin ~ 1-4 ± 0.1 
days and the SF curve is characterized by a slope /3 — 0.65 ± 0.06 
until tmin and /? — 1.37 ± 0.17 from tmin to tmax. Below t ~ 3 
days, the variability at 2.8 cm appears more pronounced than that at 
1 1 cm and 6 cm. However, Table 1 shows that the measurement un- 
certainty mo at 2.8 cm is I.l per cent, which is 3 to 4 times higher 
than at 1 1 cm or 6 cm. Although the formal significance of variabil- 
ity at 2.8 cm is higher than 99.9 per cent, the larger measurement 
uncertainty at this wavelength should lead to a cautious interpre- 
tation of the variability in this band. While in the first half of the 
2.8 cm light curve the variability is very low (m = 1.15 per cent), 
we note in the second half a continuous rising trend in flux density, 
which is the main reason for the high significance of the overall 
variability. Since a departure from this rising trend is not seen, the 
true variability amplitude could be higher and the variability time 
scale could be longer than 1.4 days. 

In summary we find that the fastest and most pronounced vari- 
ability is seen at the shortest radio wavelength (2.8 cm) with a char- 
acteristic time scale of 0.5 ± 0.05 days. Common to the three radio 
data sets (11,6, and 2. 8 cm) is also a slower variability mode, which 
appears at a characteristic timescale of 1.1 days at 6 cm and to 
1.4-1.5 days at 1 1 cm and 2.8 cm. The structure function at 2.8 cm 
differs substantially from the more similar ones seen at 6 cm and 
1 1 cm. Such different shapes toward higher frequencies were also 
seen previously in this source (Fuhrmann et al. 2008) and were in- 
terpreted as indication for the presence of physically different vari- 
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Figure 1. The broad band (radio through optical to X-ray ) variability observed in the source S5 0716+714 over the campaign period, (a) 0.2—10 KeV light 
curve of the source; (b) optical light curves of the source in B,V,R and I passbands; (c) flux density variability at 2.8 cm wavelengths; (d) the flux density 
variability at 6.0 cm from two different telescopes (Effelsberg and Urumqi) and (e) flux density variability at 1 1.0 cm wavelengths. 



ability mechanisms, as e.g. interstellar scintillation dominating at 
the longer wavelength and a with frequency increasing contribu- 
tion of source intrinsic variability. In this context, variability time 
scales which are seen in the 2.8 cm data, but are not seen in the 
1 1 cm and 6 cm data, may be related to some source intrinsic vari- 
ability and may well therefore appear also in the optical, which is 
not affected by interstellar scintillation. 

In the optical bands S5 0716+714 also displayed remarkable 
variability during the campaign period. The most dense coverage at 
optical frequencies was obtained in the V and R passbands. The SFs 
of these data are shown in Fig. [3] They show rapid variability with 
multiple cycles of rises and declines. Differences between both SFs 
are likely due to the different time-sampling in both data sets and 
the location of observing gaps within the two data trains. The B 



and I passband light curves are much sparser and so cannot reveal 
significant features in the SF analysis. We therefore focus on the 
SFs for R- and V-bands. The first break in both R and V band SF 
curves appears at a timescale of tmin ~ 0.25 ± 0.05 days. At R- 
band the SF shows pronounced maxima also at 0.50 ± 0.05 days, 
0.90 ± 0.05, 1.4 ± 0.1 days and 1.9 ± 0.1 days. At V-band the 
variations are less rapid, with a dominant time scale at 0.8 ± 0.05 
days and 1.9 ±0.1 days. The 0.5 day and the 1.4 ± 0.1 day time 
scale, seen in R-band, are not seen in V-band. 

We note that these time scales are within their measurement 
errors factor of two multiples of each other, providing a hint of 
an underlying 0.25 day quasi-periodicity (QPO). We further note 
that at 2.8 cm the shortest timescale seen in the SF is 0.5 days, 
which is very close to the second harmonic of the optical variabil- 



Multiwavelength i 



ity. This might indicate some correlation between radio and optical. 
We check the possibility for such a correlation in Section 3.3. We 
investigate the possibility for QPOs being present in the optical data 
using other techniques in Section 3.6. 

3.3 Auto- and Cross-correlations 

To quantify the correlation among the multi-frequency light curves 
of the source during the campaign period, we computed the discrete 
auto and cross correlation functions (ACF, DCF) between different 
frequencies to search for possible time lags. We followed Edelson 
& Krolik (1988) and Lehar et al. (1992) to calculate the DCF with 
details given in Rani et al. (2009). 

In Figure |4] we plot the two auto-correlation functions (ACF) 
of the 6 cm data from the Effelsberg and Urumqi telescopes. The 
differences seen in the shape of the two ACFs is explained by the 
different data-quality from the two telescopes, with a higher signal- 
to-noise ratio and smaller beam size (and therefore lower in beam 
confusion) at the larger telescope. 

In Figure [5] we show the DCF of the 6 cm data from the Ef- 
felsberg and Urumqi telescopes. 0716+714 was observed at both 
telescopes with very similar time sampling and time coverage. The 
high degree of correlation of the two data trains is obvious and con- 
vincingly demonstrate the reality of the observed intra-day variabil- 
ity. No significant time lag (rg.o/e.o = 0.0 ± 0.05) is seen between 
the two telescopes. 

Similar to the structure function analysis (see Fig.O the auto- 
correlation functions of the 1 1 cm and 6 cm Effelsberg data show 
a common shape and show similar decorrelation time-scales (Fig. 
[SJ. On the other hand, the 2.8 cm data decorrelate much faster, in- 
dicative of faster variability at this shorter wavelength (higher fre- 
quency). The ACFs of the 1 1 cm and 6 cm bands also suggest a 
regular pattern on a ~ 2 day timescale. Since this timescale is just 
half of the duration of the overall ~ 4 day time coverage of the 
radio data, we do not regard this timescale as significant. 

To search for possible timelags between the three radio bands, 
we performed a cross-correlation analysis of the Effelsberg radio 
data. Figure [7] shows the two independent cross-correlations of the 
11.0 cm and 6.0 cm light curves versus the 2.8 cm data. Formally 
we calculate a time lag of r = —(1.2 ± 0.1) day between 11.0cm 
and 2.8 cm in the sense that the 2.8 cm data lead. Such frequency 
dependence is consistent with the canonical behavior seen in AGN 
(e.g. van der Laan 1966, Marscher & Gear 1985). 

Next, we look for a possible correlation with flux variations at 
optical-radio frequencies. Figure [8] (top) may indicate a correlation 
with ~ (0.35 ± 0.15) day time lag between R and and 2.8 cm, in 
the sense that radio is leading over optical, which apparently is op- 
posite to physical expectations (the higher frequency should come 
first). At V-band the maximum correlation with 2.8 cm occurs at 
zero time lag, however, with a lower cross-correlation coefficient 
than for R-band (DCF peak at: 0.73 for R-band, 0.53 for V-band). 
Since the R- and V-bands cross-correlate very strongly with zero 
time lag (Fig.|9}, we think that the above time lag of ~ 0.35 day is 
unlikely to reflect physical reality. The R-band light curve shown 
in Fig[T]is less continuous and shows larger time gaps, while the 
V-band curve is more continuous and therefore better suited for the 
search of a possible cross-correlation. We generate flux-flux plots 
to quantify a possible correlation between optical and 2.8 cm data. 
In Figure [Tol we show the results using R-band data and in Figure 
[TT]with V-band data. Formally we obtain the following correla- 
tion coefficients (r being the linear Pearson correlation coefficient) 
and significances. From Figure [Tol R versus 2.8 cm: r =0.163 and 
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Figure 2. Variability index plotted versus frequency. At 10.5 GHz (2.8 cm) 
the open symbol is for the whole data set and the filled square represents the 
truncated data set without the flux density increase seen at the end (see text). 
The lines show a model for a thin scattering screen located at a distance of 
lOOpc (dashed line) and 200 pc (dotted line). 



58.59 per cent significance; R versus time shifted 2.8 cm: r= 0.687 
and 99.90 per cent significance. From Figure[TT] V versus 2.8 cm: 
r=0.52 and 99.92 per cent significance. However, if we remove the 
highest data point in Figure [TO] & [TT] the correlation significance 
drops below 95 per cent for R-band versus 2.8 cm while it is still 
sustained for V versus 2.8 cm. 

We therefore conclude that we cannot claim a significant de- 
tection of a correlation between optical R-band and 2.8 cm radio 
band, probably due to the limited time sampling of the R-band data. 
In the more continuous V-band, however, such a correlation may 
exist, although we regard the evidence for it as weak, and the nom- 
inal time lag as physically unlikely. 



3.4 Spectral energy distribution 

From the available multi-frequency data of this campaign we can 
construct a quasi-simultaneous SED of the BL Lac S5 0716+714. 
For the SED we averaged in time all flux density measurements 
obtained at the different observatories during the period December 
11-15, 2009. For radio through X-ray bands the length of the error 
bar reflects the strength of the variability. At GeV frequencies the 
error bars are derived from the spectral fitting. As the IDV obser- 
vations were carried out only at three radio wavelengths, we also 
included additional flux density measurements at other wavelength 
from the IRAM 30-m Pico Veleta telescope and some public re- 
sources. These measurements were obtained quasi-simultaneously 
during the period December 9-16, 2009 and will be discussed in 
more detail in a forthcoming paper on the long term variability of 
0716+714 (Ranietal. 2012). 

The SED of the source is displayed in Fig.[T2] The upper panel 
displays the v - uFi, plot, while the - Fi, plot is shown in lower 
panel. The maximum luminosity of the SED appears to be in the 
range of around 10^^ - 10^'' Hz, which, unfortunately is in a spec- 
tral region with no data. 

In the SED we also included the public available 7— ray data 
from the second AGN catalogue of the Fermi satellite (2FGL cat- 
alogue; Ackermann et al. 2011). These data result from a 2 year 
average (Aug. 2008 - Aug. 2010) of the source flux density and 
therefore are not truly simultaneous to our observing campaign. We 
notehoweverthat0716+714wasinalow7— ray state at the time of 
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Figure 3. Structure function analysis curves of S5 0716+714, with radio 
in the upper panel and optical in the lower Different symbols are used for 
different wavelengths. 




Figure 4. Auto-correlation functions of the 6.0 cm data from Effelsberg and 
Urumqi. 



our observation, so that the expected GeV flux in December 2009 
are likely be lower than the 2FGL values plotted in this SED. 



3.5 Jet Doppler factors 

Within the framework of synchrotron self-Compton models it is 
possible to constrain the Inverse Compton Doppler factor (Sic) by 
comparing the expected and observed fluxes at higher frequencies 
(Marscher 1987 and Ghisellini et al. 1993). This IC Doppler factor 
is defined as 



Sic ^ f{c()Sr. 



\ /){6-4ci) (5-3a) 
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Figure 5. Cross-correlation function of the 6.0 cm data between the Effels- 
berg and the Urumqi telescopes. 
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Figure 6. Auto-correlation functions of the 11.0 cm (circles), 6.0 cm 
(squares) and 2.8 cm (diamonds) data from Effelsberg. 




Figure 7. Cross-correlation functions of the Effelsberg data: circles denote 
11.0 cm versus 2.8 cm, and squai'es denote 6.0 cm versus 2.8 cm. 
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Figure 8. The cross-coiTelation analysis curves of optical R passband light 
curve (top) and optical V passband light curve (bottom) with all the three 
radio wavelength light curves of the source. A bin size of 0. 1 day is used in 
the correlation analysis. 



where vt is the synchrotron high frequency cut-off in GHz, Sm 
the flux density in Jy at the synchrotron turnover frequency i/m, 
Sic the observed 7-ray flux in Jy (assumed to arise from the IC 
process) at v-^ in keV, a is the spectral index of the optically thin 
part of the spectrum, 9^ the angular source size in mas and /(a) ~ 
0.14 — 0.08q. For the high energy cut-off we follow Fuhrmann 
et al. (2008) and use ut ~ 5.5 x 10^ GHz. The angular size of 
the rapidly varying region 6^ is constrained to be ^0.03 mas as 
used by Agudo et al. (2006) and consistent with size measurements 
from mm-VLBI. We adopt this value in our further calculations. 
For the radio spectrum we use an optically thin spectral index of 
Q — —0.44, and for the synchrotron turnover we adopt Sm = 8.9 
Jy and = 148 GHz (for details of spectral fitting see Rani et al. 
2012). 

The estimated limits on Sic are 18.8 (adopting Sic = 
1.77x10^^ Jy at = 1.4 keV) and Sic = 26.1 (adopting Sic 
= 1.73x10"^ Jy at = 7.2 keV). If the same calculation is ap- 
plied for the GeV fluxes and if we assume that the same pho- 
tons would be producing also the gamma-rays, then a consider- 
ably higher Doppler factor of Sic = 69 is obtained (adopting 
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Figure 9. Cross-correlation analysis curves between the optical V- and R- 
band light curves. Although a bin size of 0.01 day is used in the correlation 
analysis this curve is shown with a bin size of 0.02 day. 




Figure 10. The 2.8 cm radio flux plotted versus optical R-band flux. The 
blue symbols show the original data, the time shifted 2.8 cm radio data are 
shown in red (a 0.35 day shift was applied, see text). 



Sic = S^ = 1.95x10"" Jy at = 5.47x10^ KeVfl For a 
given Lorentz-factor 7, the maximum Doppler-factor is ~ 27. A 
Doppler-factor of up to ~ 50 therefore is consistent with the ob- 
served jet speed, which is in the range of 7 = 15-25 (Bach et al. 
2005, Fuhrmann et al. 2008, Rastorgueva et al. 2009). 

In the inverse Compton process it is usually assumed that the 
synchrotron photons from the IR-optical branch are upscattered up 



* The GeV fluxes were obtained from public available FERMI/LAT light 
curves at |http://fermi.gsfc. nasa.gov/ssc/data/access/lat/msl_ lc/| 
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Figure 11. The 2.8 cm radio flux plotted versus optical V-band flux. 



to GeV energies. If we assume that the jet is not stratified, i.e. that 
the Doppler-boosting is homogeneous over the emission region, we 
find that the observed GeV flux could also be explained via inverse 
Compton scattering from a synchrotron component with a spec- 
tral maximum in the 300-400 GHz regime. The possibility of SSC 
dominance in Compton dominated blazars has been investigated in 
detail in Zacharias & Schlickeiser (2012). 

It is also possible to obtain a limit on the Doppler factor 5 
by assuming that the high-energy 7-ray photons can collide with 
the softer radiation to produce e* pairs. The cross-section of this 
process is maximized at ~ (tt/S (see Svensson 1987 for details), 
where ctt is the Thompson scattering cross-section. This leads to a 
lower limit on 5 (following Dondi &. Ghisellini 1995): 



5 > 3.5x10 



3 (1 + zy^jl + z - ^/T + 7)V43.8t/i/.)"] 1/(4+2") 



tvar 



(2) 

where is the X-ray flux in /ijy and i^x is the corresponding X- 
ray frequency in keV, v is the GeV frequency, a is the spectral in- 
dex measured between 1 KeV and 100 MeV. Using Fx = 1.77 ^]y, 
Vx=\AA KeV, a = -1.05, we obtained 5 > 24 for 300 MeV and 
(5 12 for 1 GeV energies. Since the published GeV data are not 
simultaneous with our radio to X-ray data, the accuracy to which 
the Doppler factors can be estimated is limited. A more accurate, 
though not largely different, estimate of the Doppler boosting will 
be presented in a forthcoming paper, which is based on truly simul- 
taneous gamma-ray data (Rani et al. 2012). 

From all of the above approaches we think it reasonable to 
conclude that the blazar 0716-1-714 exhibits a high Doppler factor 
of at least 12 and most likely ^ 20 — 25, which is not typical for a 
BL Lac type object. 



3.6 Nearly periodic variations in tlie optical? 

The complete optical light curve is given in Figure [T3] with the 
different observatories labelled. It is important to note the excel- 
lent agreement between the data taken at different telescopes dur- 
ing the many periods when we had essentially simultaneous cov- 
erage. Over the course of the 7 total days of these observations 
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Figure 12. Broad band (radio to 7-ray) spectrum of S5 0716+714 over the 
campaign period. Simultaneous data are in blue color. For comparison we 
also include the non-simultaneous 7— ray data from Fermi/LAT from the 
2FGL catalogue (red symbols, see text). The top panel shows the spectral 
energy {uFu), the bottom panel flux density {Fu) plotted versus frequency. 
For the radio through X-ray bands the error bars represent here the varia- 
tions in flux density rather than the measurement uncertainty. 



the R magnitude of S5 0716+714 varied between 12.48 mag and 
13.25 mag. Strong variations were seen throughout the campaign, 
with frequent excursions of 0.2 - 0.4 mag over a few hours. Such 
changes are not uncommon for this active blazar (e.g. Montagni et 
al. 2006; Gupta et al. 2008; Rani et al. 201 1 and references therein). 
The duty cycle in previous observations was found to be unity, 
with some detectable variation every time it is observed for more 
than a few hours (Wagner & Witzel 1995, Wagner et al. 1996) and 
our results certainly agree with those earlier data. We have quanti- 
fied these variations by examining them for possible timescales and 
quasi-periodic fluctuations. 

Following Mohan et al. (2012), we employ a suite of statisti- 
cal techniques to analyze the optical light curve (Fig. II It. namely, 
Fourier periodogram, wavelet analysis, multi-harmonic analysis of 
variance (MHAoV) and Lomb-Scargle periodogram (LSP). The 
LSP and MHAoV offer a natural way of period detection when 
the data is unevenly sampled with the MHAoV also able to detect 
non-sinusoidal signals. Together, these are designed to detect pos- 
sible QPOs and test the validity of a detection using Monte-Carlo 
simulations based significance testing. This suite helps in making 
QPO detection reliable through consistent detection with signifi- 
cance testing in all techniques. In addition, the phase of existence 
and number of cycles a QPO is present during the entire length of 
the observations is obtainable from a wavelet analysis. As an exam- 
ple of this suite technique we apply it to our new data on 07 16+7 14. 

The fitting of the Fourier periodogram (top plot in Fig.ll4t fol- 
lows the procedure described in Vaughan (2005). It is well fit by a 
power law model with a slope of— 2.0±0.16). The Kolmogorov- 
Smimov goodness of fit test yields a p value of 0.99 indicating a 
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good fit. Errors on the fit parameters, the slope and normahzation 
are determined. Then, the error on the power law model is evalu- 
ated, after which model uncertainties are accounted for before plac- 
ing a 99 per cent significance contour (top plot in Fig.ll4t. There 
are no significant detections above this level. Monte-Carlo (MC) 
simulations are carried out using the Timmer & Konig (1995) algo- 
rithm to generate 6600 random light curves with a range of slopes 
within the fit parameter errors to simulate the properties of the orig- 
inal light curve. The MHAoV and LSP are determined for each of 
the generated light curves and an estimate of the number of times 
the data based MHAoV and LSP ordinates at the interesting period- 
icity are above the simulated ordinates gives a measure of the sig- 
nificance. The MHAoV (bottom plot in Fig.ll4t detects a possible 
period at 1.1 days with a probability of 94.6 per cent from the MC 
simulations based significance test. The LSP (bottom plot in Fig. 
114b detects the same period of 1 . 1 days with a probability of 80 per 
cent from the MC simulations based significance test. The wavelet 
analysis (Fig.|15t indicates a quasi-period of 0.93t{;;2| days in the 
duration between ~ 7.4 days and 9.4 days within the cone of influ- 
ence (triangular bounded region in between the two black lines in 
Fig.|15b indicating that it lasts for 2.2 cycles. By integration of the 
wavelet signal along the time axis (abscissa in Fig.llSb. the global 
wavelet power spectrum (GWPS) is obtained (bottom plot of Fig. 

Since the nominal periodicity shows a fairly large range of 
0.65 days -1.1 days in the optical data, and since it lasts only for 2 - 
3 cycles, any periodicity that may be present is certainly broad and 
weak. We can argue that a timescale of ~ 1 day has been seen in 
the optical variations during our campaign as it is supported by the 
moderate significance estimated for the MHAoV and LSP detec- 
tions. However, despite having a relatively large amount of optical 
data from telescopes at different longitudes, there were gaps during 
the campaign that were only filled in by the Swift UVOT monitor. 
These data are of lower accuracy than the bulk of the ground-based 
data and they affect the sensitivity of our tests for QPOs. 



4 DISCUSSION AND CONCLUSIONS 

We have successfully carried out a multi-wavelength campaign to 
observe the well-known blazar S5 0716+714. The optical portion 
of the data that we have reported here was obtained from 9 ground- 
based and one space-based telescopes over almost 7 days. The tar- 
get was very active, with repeated 'mini-flares' of 0.2-0.7 mag in 
amplitude and an overall variability of 0.8 mag in the optical V- 
band. Significant peak-to-peak variations of the order of 5-10 per 
cent were also seen in the radio bands but only marginal variabil- 
ity was detected in the X-rays. In the radio bands the variability 
amplitudes increase with frequency, and the variability time scales 
become shorter. We find significant correlations between the differ- 
ent radio bands as well as between the optical passbands, but no 
significant correlations between radio/optical and X-ray bands. 

The visual inspection of the optical light curve suggests the 
presence of a characteristic variability timescale. Therefore we an- 
alyzed this formally using the periodogram, LSP, MHAoV and 
wavelet techniques. A timescale in the range of 0.9 — 1.1 days is 
determined from these. Because of its moderate formal significance 
(80 per cent from LSP and 94.5 per cent from MHAoV analyses) 
and a small number of putative cycles (2.2 seen in the wavelet anal- 
ysis) a formal claim of strict (or quasi-) periodicity cannot be made. 

The structure function analysis yields for the three fastest vari- 
ability modes timescales of 0.25, 0.5 and 1.0 day in the optical 



bands. In the radio bands the three fastest variability modes ap- 
pear at 0.5, 1.1 and 1.5 days. Within the measurements' uncertain- 
ties, the variability timescales at 0.5 day and 1.0 day appear com- 
mon between the radio and optical. It may be noteworthy that these 
timescale are essentially integer multiples of the fastest optical vari- 
ability timescale of 0.25 days. It is unclear whether this is a chance 
coincidence or reflects a common physical origin. 

A comparison of the observed modulation indices at radio 
bands with the expected strength for (weak) interstellar scintillation 
(ISS) suggests that most of the observed radio IDV is in agreement 
with a ISS slab model, which is able to reproduce the observed 
variability time scale of the order of ~ 1 — 2 days. At 2.8 cm, how- 
ever, the observed variability index is larger than expected from 
the model which could be due to some 'underlying' source intrin- 
sic variability. This interpretation is supported by the detection of 
a time lag between 2.8 cm and the longer wavelengths in the sense 
that the 2.8 cm variability is leading. This is typical for AGN vari- 
ability and is commonly interpreted as a source intrinsic opacity 
effect. 

We note a possible cross-correlation between 2.8 cm and opti- 
cal V-band. Although mathematically significant, such a correlation 
must be regarded with some skepticism in view of the lack of a sim- 
ilar correlation at R-band, the unexpected sign of the putative time 
lag (radio is leading) and the limited length of the data trains. So, 
we cannot claim any significant correlation between the flux vari- 
ations at optical and radio frequencies over this 1-week IDV cam- 
paign. We further note that in their analysis of the long-term vari- 
ability of 0716-1-714 during 2007-201 1, in which 0716-1-714 shows a 
prominent flare with peak in late December 2009, Rani et al. (2012) 
find that the optical flux variations lead the radio by ~ (60 — 70) 
days. Such a time lag would exclude any direct radio-optical corre- 
lation in the 1 week long data train discussed here. 

As a matter of fact, correlated radio-optical variability of 
this source has been investigated in several earlier observing cam- 
paigns. During four weeks of multi-frequency observations in 
February 1990, Quirrenbach et al. (1991) detected a significant one- 
to-one correlation between optical/radio flux variability (see also 
Wagner et al. 1996, Qian et al. 1996). Such a correlation was not 
seen in two later radio-optical IDV campaigns, which were per- 
formed in 2000 (unpublished data) and in 2003 (Ostorero et al. 
2006, Agudo et al. 2006, Fuhrmann et al. 2008). 

To further explore the possibility of a radio-optical correlation, 
we examine the spectral characteristics of 0716-1-714 for three dif- 
ferent radio-/optical IDV campaigns. In Table[2]we summarize for 
each observing date the spectral properties of 0716-1-714, where we 
used flux densities averaged over the duration of the campaign. In 
the table we list the synchrotron peak frequency (i^max) in column 
2, the radio spectral index in column 3, and the mm-optical spec- 
tral index in column 4. We see that in 1990 the source showed a 
factor ^ 3 lower synchrotron turnover and considerably steeper ra- 
dio spectrum than in the two later campaigns, where 0716-1-714 was 
much more active and showed prominent flux density outbursts. 

It is therefore likely that the detection/non-detection of corre- 
lated radio-optical IDV relates to the activity state of the source. IN 
particular, the opacity of the emitting region in the cm-bands (e.g. 
Qian 2008) and the presence of high or low peaking flux density 
flares, which trace moving shocks in a relativistic jet at different 
separations from the jet base (e.g. Valtaoja et al. 1992), may deter- 
mine the presence of such correlations. Near and above the actual 
turnover frequency, the opacity decreases progressively, the spec- 
trum steepens and the time lag between radio and optical variations 
vanishes. The direct radio-optical correlation seen in February 1990 
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Table 2. A comparison of spectral variations over different campaign 
epochs 

Epoch i^iiiu.r <h-uiJii} Oiiiiii oiiiiciii Reference 

Feb. 1990 35±5 -().82ib0.08 -0.75±0.05 1 
Nov. 2003 90±5 -0.35±0.10 -0.88±0.04 2 
Dec. 2009 114±12 -0.44±0.10 -0.85±0.03 3 

i^max '■ synchrotron peak frequency in GHz from spectral fit 
aradio '■ Optically thin radio spectral index calculated by fitting a syn- 
chrotron self-absorbed model over a frequency range 2.7 - 230 GHz. 
otmrn-optical '■ Spectral index from 230 GHz to optical R-band. 
References : 1. Quirrenbach et al. 1991; 2. Ostorero et al. (2006) and 3. this 
paper . 

implied a small (< 1 day) time lag in the variability, which is fully 
consistent with the obsei-ved optically thin state and a time of inter- 
mittent quiescence in the overall source variability. 

The detection of QPOs with similar time scales in radio and 
optical band would restrict their physical origin to jet rather than 
to accretion disk physics. As the activity state of the source and its 
time variable opacity play a fundamental role in the search for cor- 
related radio-optical IDV, we see two ways to proceed in the near 
future: (i) observe during times of low activity, e.g.. minimal flux 
and steep spectrum, and/or (ii) observe at frequencies above the 
critical turnover frequency, which for 0716-1-714 is at ^ 100 GHz 
(in mm- and sub-mm bands). In the optically thin regime the lower 
source opacity will ensure that the observer actually looks at the 
same physical region. In order to facilitate the identification of 
QPOs, data trains have to be long enough in time so that a suffi- 
ciently large number of variability cycles can be observed. With a 
typical variability time scale of 0.5-1.5 days in 0716-1-714, at least 
5—10 duty cycles are necessary to unambiguously identify a QPO. 
This means that a continuous broad band monitoring for more than 
1 week duration is necessary. In order to separate the intrinsic vari- 
ability from unavoidable interstellar scintillation (caused by the 
small size of the emission region) multi-frequency observations are 
required, with a frequency coverage which should include the mil- 
limeter and sub-millimeter radio bands. This will help to avoid the 
effects from interstellar scintillation, which dominate at the longer 
cm-wavelengths. 



ACKNOWLEDGMENTS 

We thank the referee for constructive comments that have helped us 
to improve the paper. We acknowledge data from the 100 m Eff els- 
berg radio telescope, which is operated by the Max-Planck-Institute 
ftir Radioastronomie in Bonn (Germany). This paper made use of 
the data obtained with the Urumqi 25m radio telescope of Xinjiang 
Astronomical Observatory of the Chinese Academy of Sciences 
(CAS). The IRAM 30m Telescope is supported by INSU/CNRS 
(France), MPG (Germany), and IGN (Spain). We thank the Swift 
team for making these observations possible. The work of ACG, 
BR, RB, HG, SP, ES and AS was partially supported by Scientific 
Research Fund of the Bulgarian Ministry of Education and Sci- 
ences (Bin- 13/09 and DO 02-85) and by Indo Bulgaria bilateral 
scientific exchange project INT/Bulgaria/B-5/08 funded by DST, 
India. The work of PJW is partially supported by US NASA grant 
NNX11AB90G. The Abastumani Observatory team acknowledges 
financial support by the Georgian National Science Foundation 
through grant GNSF/ST09/521_4-320. The work of KEG was par- 
tially supported by the COST Action MP0905 "Black Holes in a 



Violent Universe" and by the Hungarian Scientific Research Fund 
(OTKA, grant no. K72515). lA acknowledges the funding support 
by the regional govenmient of Andalucia, and by the Ministerio de 
Ciencia e Innovacion of Spain through grants P09-FQM-4784 and 
AYA2010- 14844, respectively. TP acknowledges data based on ob- 
servations made with the Nordic Optical Telescope, operated on the 
island of La Pahna jointly by Denmark, Finland, Iceland, Norway, 
and Sweden, in the Spanish Observatorio del Roque de los Mucha- 
chos of the Instituto de Astrofisica de Canarias. 



REFERENCES 

Ackermann M., et al., 201 1, ApJ, 743, 171 
Agudo I., et al., 2006, A&A, 456, 1 17 

Bach U., Krichbaum T.P, Ros W, Britzen S., Tian WW., Kraus 
A., Witzel A., Zensus J.A., 2005, AA, 433, 815 
Bach U., Krichbaum T.P., Kraus A., Witzel A., Zensus J.A., 2006, 
AA, 452, 83 

BertinE., Amouts S., 1996, A&AS, 117, 393 
Bhat N. D. R., Gupta Y, Rao A. R, 1998, ApJ, 500, 262 
Beckert T., Fuhrmann L., Cimo G., Krichbaum T.P., Witzel A., 
Zensus J. A., 2002, 6th European VLBI Network Syinposium on 
New Developments in VLBI Science and Technology, Bonn, Ger- 
many, eds. E. Ros, R. W. Porcas, A. P. Lobanov, and J. A. Zensus, 
p. 79 

Breeveld A. A., et al, 2010, MNRAS, 406, 1687 

Bottcher M., et al., 2005, ApJ, 631, 169 

Bdttcher M., et al, 2009, ApJ, 694, 174 

Burrows D. N., et al., 2005, SSRv, 120, 165 

Camenzind M., & Krockenberger M., 1992, AA, 255, 59 

Chandra S., Baliyan K. S., Ganesh S., Joshi U. C, 2011, ApJ, 731, 

118 

Dondi L., Ghisellini G., 1995, MNRAS, 273, 583 

Edelson R. A., Krolik J. H., 1988, ApJ, 333, 646 

Fuhrmann L., et al. 2008, A&A, 490, 1019 

Gabanyi K. E., et al., 2007, A&A, 470, 83 

Ghisellini G., Padovani P., Celotti A., Maraschi L., 1993, ApJ, 

407, 65 

Gopal-Krishna, WiitaR J., 1992, A&A, 259, 109 

Gupta A. C, Banerjee D. R K., Ashok N. M., Joshi U. C, 2004, 

A&A, 422, 505 

Gupta A. C, Fan J. H., Bai J. M., Wagner S. J., 2008, AJ, 135, 
1384 

Gupta A. C, Srivastava A. K., Wiita P J., 2009, ApJ, 690, 216 

Hardee R .E., Walker R. C, Gomez J. L., 2005, ApJ, 620, 646 

Heidt J., Wagner S. J., 1996, A&A, 305, 42 

Kalberla P. M. W., Burton W. B., Hartman D., Amal E. M., Bajaja 

E., Morras R., Poppel W. G. L. 2005, A&A, 440, 775 

Kato S., Fukue J., 1980, PASJ, 32, 377 

Kolesnikova D. M., Sat L. A., Sokolovsky K. V., Antipin S. V., 
Samus N. N., 2008, AcA, 58, 279 

Kolesnikova D. M., Sat L. A., Sokolovsky K. V., Antipin S. V., 

Belinskii A. A., Samus' N. N., 2010, ARep, 54, 1000 
Krichbaum T. P., Kraus A., Fuhrmann L., Cimo G., Witzel A., 
2002, PASA, 19. 14 

Larionov V. M., et al., 2008, A&A, 492, 389 

Lehar J., Hewitt J. N., Burke B. F., Roberts D. H., 1992, ApJ 384, 

453 

Liu X., Song H.-G., MarchiU N., Liu B.-R., Liu J., Krichbaum 
T. R, Fuhrmann L., Zensus J. A., 2012, A&A, 543, A78 
Mangalam A. V., Wiita R J., 1993, ApJ, 406, 420 



Multiwavelength Intraday Variability of the BL Lac S5 0716+714 13 



Marchili N., et al., 2010, A&A, 509, A47 

Marscher A. P., 1987, in Zensus J. A., Pearson T. J., eds, Superlu- 

minal Radio Sources, Cambidge University Press, p. 280 

Marscher A. R, et al., 2008, Natur, 452, 966 

Marscher A. R, Gear W. K., 1985, ApJ, 298, 1 14 

Marscher A. R, Gear W. K., Travis J. R, 1992, in Valtaoja E., Val- 

tonen M., eds. Variability of Blazars, Cambidge University Press, 

p. 85 

Mohan, R, Mangalam., A., Gupta, A. C. & Wiita, R J., 2012 (in 
preparation) 

Montagni R, MaselU A., Massaro E., Nesci R., Sclavi S., Maesano 
M.,2006, A&A, 451,435 

Nilsson K., Pursimo T., Sillanpaa A., Takalo L. O., Lindfors E., 
2008, A&A, 487, L29 

Nowak M. A., Wagoner R. V., 1992, ApJ, 393, 697 

Ostorero L., et al., 2006, A&A, 451, 797 

Pollock J. T., Webb J. R., Azarnia G., 2007, AJ, 133, 487 

Poon H., Fan J. H., Fu J. N., 2009, ApJS, 185, 511 

Poole T. S., et al., 2008, MNRAS, 383, 627 

Qian S. J., Quirrenbach A., Witzel A., Krichbaum T. P., Hummel 

C. A., Zensus J. A., 1991, A&A, 241, 15 

Qian S.-J., Li X.-C., Wegner R., Witzel A., Krichbaum T. R, 1996, 
ChA&A, 20, 15 

Qian S.-J., Witzel A., Krichbaum T. R, Wagner S. J., 1995, 

ChA&A, 19, 522 

Qian S.-J., 2008, ChJAA, 8, 219 

Qmrrenbach A., et al., 1991, ApJ, 372, L71 

Raiteri C. M., et al., 2003, A&A, 402, 151 

Raiteri C. M., et al., 2005, A&A, 438, 39 

Raiteri C. M., et al., 2007, A&A, 473, 819 

Raiteri C. M., et al., 2008a, A&A, 480, 339 

Raiteri C. M., et al., 2008b, A&A, 485, L17 

Raiteri C. M., et al., 2008c, A&A, 491, 755 

Rani B., Wiita R J., Gupta A. C, 2009, ApJ, 696, 2170 

Rani B., et al., 2010a, MNRAS, 404, 1992 

Rani B., Gupta A. C, Joshi U. C, Ganesh S., Wiita R J., 2010b, 

ApJ, 719, L153 

Rani B., Gupta A. C, Joshi U. C, Ganesh S., Wiita R J., 2011, 

MNRAS, 413, 2157 

Rani B., et al., 2012 (in prepration) 

Rastorgueva E.A., Wiik K., Savolainen T., Takalo L.O., Valtaoja 

E., Vetukhnovskaya Y.N., Sokolovsky, K.V., 2009, AA, 494, 5 

Rickett B. J., 1990, ARA&A, 28, 561 

Roming R W. A., et al., 2005, SSRv, 120, 95 

Sagar R., Gopal-Krishna, Mohan V., Pandey A. K., Bhatt B. C, 

Wagner S. J., 1999, A&AS, 134, 453 

Simonetti J. H., Cordes J. M., Heeschen D. S., 1985, ApJ, 296, 46 
Sokolovsky K. V., 2009, Peremennye Zvezdy Prilozhenie, 9, 9 
Sokolovsky K., Lebedev A., 2005, in Simon A., Golovin A., 
eds, 12th Young Scientists' Conference on Astronomy and Space 
Physics, Kyiv Univ. Press, Kyiv, p. 79 

StaUn C. S., Gopal-Krishna, Sagar R., Wiita R J., Mohan V., 

Pandey A. K., 2006, MNRAS, 366, 1337 

Stetson R B., 1987, PASR 99, 191 

Svensson R., 1987, MNRAS, 227, 403 

Timmer, J., & Konig, M. 1995, A&A, 300, 707 

Urry C. M., Padovani R, 1995, PASR 107, 803 

Valtaoja E., Terasranta H., Urpo S., Nesterov N. S., Lainela M., 

Valtonen M., 1992, A&A, 254, 71 

van der Laan H., 1966, Nat, 211, 1131 

Vaughan S., 2005, A&A, 431, 391 



Villata M., Raiteri C. M., Lanteri L., Sobrito G., Cavallone M., 

1998, A&AS, 130, 305 

Villata M., et al., 2000, A&A, 363, 108 

Villata M., et al., 2006, A&A, 453, 817 

Villata M., et al., 2008, A&A, 481, L79 

Villata M., et al., 2009a, A&A, 501, 455 

Villata M., et al., 2009b, A&A, 504, L9 

Wagner S. J., 1992, in Valtaoja E., Valtonen M., eds. Variability 
of Blazars, Cambidge University Press, p. 346 
Walker M. A., 1998, MNRAS, 294, 307 
Wagner S. J., Witzel A., 1995, AliA&A, 33, 163 
Wagner S. J., et al., 1996, AJ, 111, 2187 

Wu J., Peng B., Zhou X., Ma J., Jiang Z., Chen J., 2005, AJ, 129, 
1818 

Wu J., Zhou X., Ma J., Wu Z., Jiang Z., Chen J., 2007, AJ, 133, 
1599 

Zacharias M., Schlickeiser R., 2012, MNRAS, 420, 84 

Zhang X.-H., Bao G., 1991, A&A, 246, 21 

Zijlstra A. A., van Hoof R A. M., Periey R. A., 2008, ApJ, 681, 

1296 



14 Alok C. Gupta et al. 




Figure 13. R passband light curve of the blazar S5 0716+714. Different symbols and colors denote for data from different observatories. Color corrected 
V-band data from Swift were added to complement the time coverage. 
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Figure 14. Top plot: periodogram on a log-log scale. The yellow solid line indicates the white noise level which is the contribution of photon shot noise at all 
frequencies. The green solid line indicates the best fit power law to the log-log periodogram. The dashed lines above and below the best fit line indicate the 
errors on the power law model. The upper yellow dashed line indicates the 99 per cent significance contour, taking into account errors on the power law as 
well as model uncertainties. Bottom plot: Combined plot showing the LSP, MHAoV and the GWPS on a log-log scale with a black veilical line at the position 
of the periodicity of 1.1 day. The MHAoV and LSP indicate a feature at 1.1 days with significances of 94.6 % and 80 % respectively. The GWPS shows a 
minor feature at 0.93 days just to the right of this. All the curves are normalized such that the total power summed over all frequencies between the sampling 
and the Nyquist frequency are equalized. 
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Figure 15. Wavelet contour plot beside which is an intensity scale indicating the color scaling of the wavelet power. Contour plot abscissa: time duration of 
observation; ordinate: periodicities searched for in the detection of quasi-periodic signal. A broad feature between JD 2455170 + 7.4 and 9.4 days is indicated 
inside a contour drawn between periodicities of 0.65 days and 1.08 days (horizontal red lines) within the cone of influence. This contour indicates the wavelet 
power above 90% confidence in this region. A periodicity of 0.93to2| is indicated which lasts for 2.2 cycles. 



